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Asymptotic behavior of the distribution functions of local 
quantities in disordered conductors is studied in the weak dis- 
order limit by means of an optimal fluctuation method. It is 
argued that this method is more appropriate for the study of 
seldom occurring events than the approaches based on non- 
linear cr-models because it is capable of correctly handling 
fluctuations of the random potential with large amplitude as 
well as the short-scale structure of the corresponding solutions 
of the Schrodinger equation. For two- and three-dimensional 
conductors new asymptotics of the distribution functions are 
obtained which in some cases differ significantly from previ- 
ously established results. 



I. INTRODUCTION. 

It has been well known for at least a decade [EJ that an 
adequate description of the Anderson localization transi- 
tion in disordered mesoscopic conductors must necessar- 
ily be formulated in terms of the full distribution func- 
tions of conductance g and/or other quantities charac- 
terizing the sample. As a consequence, even well into 
the metallic regime it should be possible to observe the 
onset of localization by studying deviations of asymp- 
totic tails of the distribution functions from their be- 
havior in the infinite conductance limit. Such a study 
was first performed in Ref. B within the framework of 
the replica-based nonlinear a- model ||. It was demon- 
strated that the tails of the distributions of such quan- 
tities as conductance, local density of states, current 
relaxation times, etc. are all described in two spa- 
tial dimensions by rather similar logarithmically normal 
asymptotes whereas Gaussian distributions would be ex- 
pected for g — > oo. Recently, it was discovered Q that 
these asymptotes can be obtained in a more straightfor- 
ward and elegant manner using the optimal fluctuation 
method in conjunction with the supersymmetric variant 
of the nonlinear a- model |5j. In a series of subsequent 
publications the method was employed to study statis- 
tics of eigenfunction amplitudes in weakly localized two- 
dimensional conductors Q and then extended to systems 
in three spatial dimensions |?] ^] . 

The common feature of all the results obtained so far 
with the use of various cr-models is that asymptotic be- 
havior of a distribution function V (t) for large t has the 
following form for two- and three-dimensional systems: 



Here d is the dimensionality of the system and t may 
have the meaning, e.g., of the current relaxation time, or 
the normalized local density of states at the Fermi energy 
p (Ep, r) /vd, or the normalized square of the local wave 
function amplitude V \ ipE ( r )| 2 of an eigenstate with en- 
ergy E. Vd is the average d-dimensional density of states 
and V is the volume of the sample. One notable excep- 
tion is the e~ ln law obtained in Ref. for the distri- 
bution of wavefunction amplitudes in three-dimensional 
samples. We will discuss a possible origin of this differ- 
ence in Section IV. 



The coefficients Cd in general depend on the strength 
of disorder and, for d = 2, on the sample size L. In the 
two-dimensional case the calculations based on cr-models 

give C2 = — ^— [§j6|-||], where pp is the Fermi mo- 
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mentum, D = — Ivp is the electron diffusion constant, I is 

the mean free path (which is assumed to be much larger 
than the electron wavelength pp 1 ), and vp is the Fermi 
velocity. (3 is a numerical coefficient which, depending on 
the symmetry of the ensemble of random potentials takes 
the following values: /3 = 1 in time-reversal-invariant sys- 
tems (the orthogonal symmetry class), and (3 = 2 when 
symmetry with respect to time reversal is completely bro- 
ken (the unitary symmetry class). In Ref. H only the 
unitary ensemble was considered. We will argue below 
that this value of C2 can only be regarded as an order- 
of-magnitude estimate. 

In the one-dimensional case all eigenstates are local- 
ized, and the distribution of wave function amplitudes 
has a simple Gaussian form (see below, Eq. (|3l|)). Cur- 
rent relaxation times t, on the other hand, are character- 
ized by a much broader logarithmically normal distribu- 
tion 



V (t) ~ exp (-C*i In 2 t) , 



(2) 



V{t) ~ exp(-C d ln d t 



(1) 



where C\ = 1/2L for a sample of length 2L ]g| JTT|] . 

In Refs. [0 and || saddle-point solutions of the su- 
persymmetric nonlinear cr-models were obtained for the 
three-dimensional case. However, conventional nonlinear 
cr-models used in Refs. [^|4|^-|8| are low- "energy" effec- 
tive field theories in which the role of energy is played 
by the diffusion operator Z?V 2 . As such they are only 
applicable to the description of phenomena that can be 
characterized as diffusive - the mean free path I is the 
smallest length scale in these theories. Formally this is 
expressed as a requirement that the scale of spatial vari- 



1 



ations of the fundamental variables of the theory - the 
Q-matrices - must be much larger than I. The optimal 
fluctuations of the Q-matrices computed in Refs. [0J|] 
were found to vary rapidly over distances ~ I in three- 
dimensional systems, which made it impossible to obtain 
rigorous results. The coefficient C3 was estimated in Ref. 
H to be of the order of [ppl) 2 ■ A similar calculation in 
Ref. § led to the same estimate for C3 but a log-normal 
rather than e~ C3 ln * functional form of the distribution 
function. 

In an attempt to overcome the limitations of diffusive 
cr-models and account exactly for the spatial variation of 
the optimal fluctuation in 3D systems at the scale of the 
mean free path, a generalized version of the model (ballis- 
tic nonlinear cr-model) was introduced by Muzykantskii 
and Khmelnitskii in Ref. and used in Ref. || to com- 
pute the distribution of current relaxation times. The 
functional form of Eq. (TP) was reproduced in that calcu- 

lation and a numerical value — = (vpI) of Ca was found. 

9^3 

We will argue based on the calculations presented be- 
low that all the results for the three-dimensional systems 
quoted above contain serious errors. 

In this article we propose an alternative method for the 
investigation of distribution function asymptotics. In- 
stead of integrating out the disorder degrees of freedom 
and then looking for a saddle point of the resulting effec- 
tive field theory (the nonlinear tr-model), we suggest that 
large-t behavior of the distribution functions is governed 
by a saddle point of the original theory in which the dis- 
tribution functions are expressed as functional integrals 
over both electronic and disorder degrees of freedom. In 
other words, such a saddle point corresponds to an op- 
timal fluctuation, i.e. the highest-probability configura- 
tion of disorder that lets the electronic eigenstates at a 
given energy E have the desired property - for example, 
an anomalously large intensity V \ tpE ( r )| 2 at some point 
r. In our view, this approach possesses the advantage 
of being applicable to systems of arbitrary dimensional- 
ity as well as to disorder models more general than the 
commonly used model with a Gaussian distribution of 
the potential fluctuations. It is also conceptually simpler 
and makes the physical origin of the results much more 
transparent. In many respects our approach is similar to 
the ideas utilized in Refs. [|l2|,^3). We were able to repro- 
duce the general ln V ~ — ln d t behavior of the distribu- 
tion functions for d = 2, 3 as well as the log-normal or 
Gaussian (depending on the precise definition of t) form 
of V for d = 1. In two- and three-dimensional cases we 
obtain new values for the coefficients Cd correcting what 
we believe to be the errors in their previously published 
estimations. 

The rest of the paper is organized as follows. In Section 
H we briefly describe the method and present the main 
results. A detailed derivation of the basic equations and 
the analysis of the saddle point solutions is deferred until 
Section III. In Section IV we will compare our results 



with those obtained by employing nonlinear cr-models. 
A short summary and a list of open questions can be 
found in Section |y|. 



II. ASYMPTOTES OF THE DISTRIBUTION OF 
WAVE FUNCTION AMPLITUDES. 

A typical wave function in a metallic sample is spread 
more or less uniformly throughout the sample volume 
V ~ L d so that its amplitude does not differ much from 
the average value of l/W anywhere in the sample. In 
the metallic (or weakly localized) regime such states ac- 
count for the bulk of the distribution of the local den- 
sity of states p (r) or of the current relaxation times t. 
The chances of observing anomalously large values of 
these quantities are related to the probability of finding 
an "anomalously localized state" . Such a state would 
be characterized by an amplitude reaching a value much 
larger than the average at some point r inside the sample. 
In what follows we will concentrate on the distribution of 
wave function amplitudes. Other distributions, such as 
those of current relaxation times, local density of states, 
etc. can be more or less straightforwardly derived within 
the framework of the same formalism. 

The problem can be formulated in the following way. 
Let us consider a spherical (in 3D) or a disk-shaped (in 
2D) conductor of radius L and compute the probability 
that an eigenstate ip of energy E (which we take to be 
equal to the Fermi energy Ep) has an amplitude ytfV 
in the center of the sample (r = 0) , with i > 1 |l4j] . The 
distribution of the disorder potential is assumed to be 
Gaussian, 



W [U (r)] = Afu exp J U 2 (r) dr 



(3) 



where Mjj is the normalization constant and r is the mean 
free time. Such a probability is then naturally expressed 
as 



p{t) = (s (V|V(o)| 5 



(4) 



where (. . .)u denotes the averaging with the weight W 
over all possible configurations of U, and ip (r) is the so- 
lution of the Schrodinger equation with the Hamiltonian 
H = H a + U (r) and energy E. Hq = p 2 /(2m) is a 
Hamiltonian of free particles with a mass m, and p is the 
canonical momentum operator. 

The requirement that E must be an eigenvalue of H 
was enforced in Ref. jf| explicitly by introducing the cor- 
responding S- function into the definition of V . In the 
alternative approach proposed here - the direct optimal 
fluctuation method - this requirement is easier to im- 
pose, if necessary, at a later stage in the calculations 
through appropriate boundary conditions for the saddle 
point equations flq]. 
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Rewriting Eq. ^ as a constrained functional integral 
and introducing Lag ran ge multipliers to enforce the con- 
straints (see Section III for details) one can demonstrate 
that the resulting "action" A possesses a saddle point, 
and the leading contribution to ln'P (t) is given by — ■> 
the value of the "action" at that point. With exponential 
accuracy the results are: 



V (t) ~ exp (-k (p F l) In 3 t) , (3D) , 



V (t) ~ exp 



(5a) 



(5b) 



where n is a number which is approximately equal to 
3 x 10~ 3 . The distance ro in Eq. ( J5b| ) is of the order of 
the electron wavelength p F . 

In both two- and three-dimensional cases, our answers 
have a different /-dependence compared to the results 
obtained in Refs. In the three-dimensional case 

ln'P is proportional to the first power of pf I while all 
the results derived from the nonlinear cr-model in three 
dimensions lead to a (j>fI) 2 dependence. Similarly, in 
the two-dimensional case the logarithm in the denom- 
inator contains L/tq instead of the L/l established in 
Ref. Q and reproduced with minor modifications in Refs. 
[^-^|. While in the two-dimensional case this difference 
leads only to small corrections, the answer for a three- 
dimensional sample indicates a substantially increased 
probability to observe an anomalously large eigenstate 
amplitude as compared to the previous results. 

These asymptotes correspond to the optimal configu- 
rations of the potential U (r) which are essentially Bragg 
mirrors, 



U(t) cx 



sin 2ppr 



r d-l 



(6) 



as r — > co. The r — > behavior is more complicated. In 
one- and two-dimensional systems the r — > region does 
not contribute to ln'P in the leading in lnt order. In 
three dimensions the shape of the optimal fluctuation at 
small distances is a narrow potential well surrounded by a 
wider potential barrier (Fig. such that the whole struc- 
ture supports a narrow resonance in the s-wave channel. 
Anomalously large values of t are achieved by combin- 
ing this local resonance with Bragg reflection at larger 
distances. 

The role of local resonances in producing large val- 
ues of t for any d in the interval 2 < d < 4 can be 
illustrated by the following simplified version of the di- 
rect optimal fluctuation method. Let us approximate the 
shape of the resonance-producing potential as U (r) = 
UqO (r — a) 9 (b — r), where 9 (r) is the step function, and 
Uq (> Ep) and b are optimization parameters. The con- 
stant a is determined from the condition that there is 
only one resonance in the potential well and it has a 
given energy Ep (a ~ pp 1 = (2mEp) 1 ^ 2 ). The wave 



function of the resonant state in such a potential is char- 
acterized by exponential decay in the interval a < r < b. 
Assuming that the amplitude of the wave function for 
r > b is close to its average value 1 / \/V we can estimate 

t = Vip 2 (0) - exp ^2^2m (U Q - E F ) (b-a)\. Antici- 
pating the final result b ~ Pp 1 \nt we can neglect a in the 
exponent. We then have 



Uq « Ep 



1 



2m V 2b 



lnt 



and for a<d 



A = 



TTV d T 



J U 2 (r) dr 



TTV d T 



2 U d 



V d U'b 



(7) 



(8) 



where Vd 



T d/2 



is the volume of a d-dimensional 



r W 2 + 1 ) 

sphere of unit radius. Substituting Eq. (|7|) into Eq. 
and minimizing with respect to b we find 



-1 / 4 - d , 



(9) 



and it follows then that Uq = ^Epj (4 — d). The resulting 
estimate for the asymptotic behavior of the distribution 
function reads 



■(t) ~ exp {- K d (p F l) ln d t} , 



1 ( 1 



2d V 4tt J V 4 - d 



v d 2 



(10a) 



(10b) 



In the three-dimensional case the coefficient in front of 

n 1 

ln J i is (p F l) fa 0.032 (p F l). It differs only by a 

number from a more rigorous estimate quoted in Eq. ( pq ) 
suggesting that local resonances play an important role 
in the three-dimensional case. As it should be, the coeffi- 
cient K3 0.032 is larger than k since it would be naive 
to expect a box-shaped potential to be optimal. 

In the two-dimensional case, large values of t are not 
optimally produced by local resonances: according to 
Eq. (pal) the coefficient in front of In 2 t is proportional 



to 



i.e. it depends on the sample size. 



The hypothesis that log-normal (in one and two di- 
mensions) and e~ Ca ln * (in 3D) asymptotics of the dis- 
tribution functions may be produced by Bragg reflection 
locking the states in (or out of) certain regions of the sam- 
ple was put forward by Muzykantskii and Khmclnitskii 
in Ref. || . They have explicitly demonstrated that a po- 
tential of the form U (x) = (pp ln (fA) /mL) cos (2p F x), 

where A = is the mean level spacing, traps an elec- 

VdV 

tron of energy Ep in a one-dimensional sample of length 
2L for time t. The probability of finding such a configu- 
ration is 



3 



r(t) 



exp 



2L 



In 2 (tA) 



(11) 



which is the same as the log-normal long-time asynrp- 
totics of the exact solution for the conductance 
Our results confirm this hypothesis in higher-dimensional 
samples. The most striking feature of our results is the 
fact that a single configuration of the potential - i.e. a 
single compact region of the configurations space - seems 
to be responsible for the large-t tails of the distribution 
functions. The alternative would be for a large number of 
more or less equal- weight configurations to be able to pro- 



duce the required wave function amplitude y — , and the 

total probability V (£) would then be given by the sum of 
the weights of these configurations. A strong argument 
against such an alternative is the fact that the results 
obtained by simply singling out the most probable type 
of configurations can lead to higher probabilities (signif- 
icantly higher in the 3D case) than the cr-model results 
which correspond to a sum over all configurations |l7j . 



III. THE DIRECT OPTIMAL FLUCTUATION 
METHOD. 

In order to enforce the condition that i\) in Eq. (Q) is a 
solution of the Schrodinger equation we rewrite Eq. (^) 
formally as 



V (t) = J VU (r) W [U (r)] A^ [U] J V^ (r) 

x l[5((H -E)^ (r'))6 (V ^ (0)| 2 - t) , (12) 



where M^, [U] is a [/-dependent normalization constant 
defined by 



N^U] J ViP (r) -£)</>(r') 



= 1. 



(13) 



Introducing auxiliary variables \ (r) and A associated re- 
spectively with the first and the second ^-functions in Eq. 
( |l2| ) enables us to represent the distribution function as 
an unconstrained functional integral: 



V (t) =N V VU (r) / V^ (r) / V 



dX 
2T 



X(r) 
2ir 



(14a) 



A [U, V>, X , A] = 2^1 J dv U 2 (r) - zA (y \i> (0)| 2 - t) 
-i J dr X (v)(^- + U (r) -E\iI> (r) . (14b) 

The n ext step is to make a saddle-point approximation 
in Eq. (14a) utilizing the fact that t is assumed to be 



large. The stationarity condition for the action Eq. ( |l4b| ) 
leads to the following set of equations: 



ttvcltU (r) - ix (r) ip (r) = 0, 



2m 



+ U (r) - E )ip (r) = 



(15a) 



(15b) 



^— + U (r) - E ) X (r) + 2A VVtS (r) = 0, (15c) 



v\ip(o)r -t = o. 



(15d) 



In deriving the saddle-point equations we neglect the 
[/-dependence of the normalization constant A/^,. The 
reasoning leading to this approximation is explained in 
detail in the Appendix. 

Below we will investigate only those solutions of Eqs. 
(U) which possess full rotational symmetry. While it 
is likely that even for rotationally invariant boundary 
conditions there exist solutions which break rotational 
symmetry we will assume that they are not optimal, i.e. 
that they correspond to extremum points other than the 
global minimum of A. The mean free time enters these 
equations only through the 'kv^t factor in Eq. (15a). It 
can be absorbed into redefinitions of x ( r ) an( l A an( l 
therefore the above assumption does not depend on the 
disorder strength. 

Eliminating X in favor of U and ip & n d switching to di- 
mensionless variables v (r) = U Q r ^ , r — pp \r\, we arrive 
(in the absence of magnetic field) at the following system 
of equations (^ (.) = (.)'): 



y (r) = In' (r {d - 1)/2 ^ (r)) 



y'(r) + y 2 (r) + ^S d , 2 + 1 = v (r) 



2y (r) (r^v (r)) - (r^v (r))' = A, 



where 



2i\t 



Pf 



-KV d T ElSd 



(16a) 



(16b) 



(16c) 



(17) 



5(1,2 is the Kronecker symbol and Sd — dVd is the ar ea o f 
a c?-dimensional sphere of unit radius. Equations (16a) 
and (16b) together are equivalent to the Schrodinger 
equation for the radial part of the wave function, while 
Eq. ( |16c| ) provides the nontrivial connection between the 
wave function and the potential that is necessary to 
achieve a given value of t at a minimal "cost" in terms of 
the weight function W. A can be viewed as a parameter 
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effecting a Legendre transformation of A, so (A, t) - al- 
ternatively, (A,lni) - form a pair of conjugate variables. 



It follows from Eq. (16c) that A must be real, so the 
saddle-point value of A is purely imaginary. The saddle- 
point action Ao is expressed in terms of the solution v (r) 
of Eqs. © as 



-KVT E% 



„d-l„,2 



(r) dr. 



(18) 



The choice of boundary conditions to Eqs. ([16]) is 
not entirely straightforward, and we would like to mo- 
tivate it by appealing to an analogy with ordinary one- 
dimensional Schrodinger equation. In the variational for- 
mulation of quantum mechanics the Schrodinger equation 
appears as a stationarity condition of a functional 



T : 



dxip' 



(*)( 



H-E) ip(x), 



(19) 



where H is the Hamiltonian, E is the particle en- 
ergy, and ip must satisfy two boundary conditions (e.g., 
4> (0) = tp (L) — 0) and also the normalization condi- 
tion f Q L dx\ip\ 2 = 1. The overall number of conditions for 
this second-order equation is three, which makes it over- 
defined. As a result, solutions exist only for a discrete set 
of values of E. One of these values would correspond to 
the ground state and thus to an absolute minimum of J- ', 
while the rest would simply correspond to its stationary 
points. In the L — > oo limit the set of allowed E values 
becomes infinitely dense so that E effectively becomes a 
continuous variable. 

An alternative way to look at the Schrodinger equation 
(often useful in numerical calculations |l^]) is to consider 
E as an extra unknown function satisfying a trivial equa- 
tion 



dx 



(20) 



One then has effectively three first-order equations with 
three extra conditions imposed on them. However, now 
these equations are nonlinear, and the three conditions 
do not specify the solution uniquely. Rather, there ex- 
ists a family of solutions corresponding to various quan- 
tum states of the system, one of which minimizes the 
functional T . Again, in the L — > oo limit the family of 
solutions becomes continuous. 

Returning to the problem of specifyin g a s et of condi- 
tions for Eqs. ( |16|) w e see that equation (16c) is a gener- 
A plays the role of an eigenvalue, 
dH) with A' = 0, 
The follow- 



alization of Eq. ( |20[ ) 
so we can formally supplement Eqs. 
increasing the number of equations to four, 
ing set of conditions must be imposed on this system of 
equations. 



First, a boundary condition is provided by Eq. (15d) 



(21a) 



The normalization requirement 

Sd [ r^drifj 2 (r) = 1 
J o 



(21b) 



provides the second condition. 

Another boundary condition is derived from regularity 
requirements on v (r) . Since v (r) can have at most a 
square integrable singularity as r — > 0, it follows that in 
three dimensions 



lim (VW>' M) = 0. 



(21c) 



In two dimensions we simply have a condition that ip' 
must be finite at r — 0, and for a one-dimensional system 
the corresponding requirement is ip 1 (0) = 0. Finally, a 
fourth condition may be provided by another boundary 
condition on ip, e.g. 



^ (L) = 0. 



(21d) 



However the problem of minimizing A is not solved 
uniquely by the set of equations ( |l6| ) together with the 
conditions (^l|). As in the ordinary quantum mechan- 
ical case, we are faced with a family of solutions, only 
one of which actually corresponds to the absolute min- 
imum Ao of the action A. Finding this "ground state" 
solution requires a parametrization of the whole discrete 
family of solutions which presents a difficult task. It 
can nevertheless be somewhat simplified by employing 
the following observation. For large enough L the set of 
solutions of Eqs. ( |l6| ) may be approximately considered 
conti nuous. It is then possible to replace the condition 
( [2 lc| ) with a free boundary condition. The result will 
be that the family of allowed solutions will now admit a 
continuous parametrization, simplifying the task of de- 
termining the minimum of the functional A. In practice 
this parametrization depends on the dimensionality of 
the sample, and we will consider the cases of d = 1, 2 
and 3 separately. Employing a continuous parametriza- 
tion which effective ly co rresponds to relaxing the bound- 
ary condition Eq. ( |21d| ) means not enforcing the condi- 
tion that Ep is an eigenvalue of the Hamiltonian. As 
was already mentioned p3] such an approximation leads 
only to errors that are beyond the exponential accuracy 
with which asymptotes of V (t) are determined in the 
optimal fluctuation method and which can therefore be 
safely neglected. 



A. One-dimensional wire. 

As a "toy" model, we consider first a purely one- 
dimensional disordered wire of length 2L. All the eigen- 
states in this case are exponentially localized with a lo- 
calization length being of the order of the mean free path 
I. A typical amplitude of such a state in its domain of 
localization is l/v2£, so the optimal fluctua- 

tion method cannot be applied to the computation of 
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the distribution function V (t) for t < L/l, Amplitudes 
larger than that are due to the states that are character- 
ized by an anomalously small localization length and the 
probability of finding such states is determined by the 
probability of the corresponding optimal fluctuation. 
The system of Eqs. (fL6|) takes the form 



v (x) = 1 + y 2 (x) + y 1 (x) , 

2y (x) v (x) — v ' (x) = -XSgnx, 



(22) 



where x is the dimensionless coordinate along the wire, 
and primes denote differentiation with respect to x. Note 
that when d = 1, y is simply the logarithmic derivative 
of the wave function. Eliminating v (x) and integrating 
once, we arrive (for x > 0) at 



(y'f = (y 2 + l) +Xy + C. 



(23) 



Here C is a constant of integration which parametrizes 
the stationary points of A as described at the end of the 
preceding subsection. Its value is fixed by a minimization 
procedure which is formally expressed as 



min„4(C) 

c 



(24) 



Although Eq. ( p3| ) is exactly integrable in terms of el- 
liptic functions, it is sufficient to make a perturbative 
expansion in A and C in order to recover the leading or- 
der dependence on t. Rewriting Eq. (|2^) as 



dx 



dy 

R(yY 



R{y) = ^(y 2 + lf + Xy + C, (25) 



we obtain the period T of the wave function oscillations 

dy 



T 



and the logarithm of the ratio of the wave function am- 
plitudes at x — and x = L: 



4> (L) T 



ydy 
R{y)' 



(27) 



We assume that the length of the sample corresponds 
to an integer number of periods, so that if y (0) = 
(from symmetry arguments) it follows also that y (-L) = 
y (L) = 0. The relative error introduced into Eq. (^) by 
neglecting the fractional part ofL/TisO(l/i). 

It is convenient to use the logarithm in the l.h.s. of 
Eq. @, which we will denote as ln#, as a parameter of 
the distribution function. It can be easily seen that the 
contribution to V of the Jacobian of the transformation 
from 9 to t is only a prefactor which does not change 
the leading order term in the exponent. Expanding the 
r.h.s. of Eqs. ( p6| ) and ( |27| ) in both A and C we obtain 
-ttA/16 = (T/2L)ln0 and T — 2-k + 0(C). For the 
action A we then find to the lowest order (with Sd = 2) 



El 2irL 

A = TIV X T — - 

Pf 32i 



(A 2 + 3C 2 ) 



(28) 



It is thus obvious that, at least to the lowest order in 1/L, 
the minimum of A is achieved when C = 0. Combining 
the results together we have for the asymptotic tail of the 
distribution function with exponential accuracy 



V (0) ~ exp 



2 4m 2 < 



(29) 



From Eq. ([25|) and the fact that the change in the wave 
function amplitude over one period (increase for x < 
or decrease for x > 0) does not depend on x, we deduce 
the exponential form of the wave function envelope 



ip (x) cx e 



• ln< 



(30) 



Normalizing the wave function to 1 and using V = 2L 
we obtain t — 2Lip 2 (0) = 2 (ln6>). We thus finally arrive 
at 



V (t) - exp 



I 

2L l 



(31) 



As explained at the beginning of this subsection the ap- 
plicability of the above formula is restricted to the re- 
gion L/l < t < ppL (the second inequality ensures the 
validity of the expansion in A and C). As required, the 
corresponding value of the localization length L/ In 8 is 
much smaller than I. 

This model calculation can be used to illustrate the 
following points. First, to the lowest order in A we have 



y (x) — cot (x + ip) + — sin 2 (x + tp) . 



(26) where ip is a phase shift, and 



v (x) = — A sin 2(x + ip) , 



(32) 



(33) 



so the shape of the optimal configuration of the potential 
is indeed a Bragg mirror. 

Second, up to a constant, 9 2 / A can be identified with t 
the delay time in the electric response of an open sam- 
ple at Ep ||, and the distribution in Eq. ( |29| ) is seen to 
coincide with the exact answer of Altshuler and Prigodin 
|ll"| quoted at the end of Section |J. Third, the sign 
of A is opposite to that of In 6, so negative values of A 
correspond to the wave function amplitude decreasing 
from the center of the sample outwards. The same will 
hold true in two and three dimensions. The amplitude of 
the oscillating potential is constant throughout the sam- 
ple. A similar observation was made in Ref. 0] — in a 
one-dimensional wire the gradients of the supersymmet- 
ric density matrix corresponding to the optimal solution 
were found to be x-independent. 



G 



B. Two-dimensional conductors 

In dimensions higher than one the system ( |l6| ) cannot 
be integrated exactly in terms of the standard elemen- 
tary or special functions. Nevertheless, numerical meth- 
ods combined with the asymptotic analysis allow us to 
fully investigate the behavior of its solutions in various 
regimes. In 2 dimensions the solution can be obtained 
as an asymptotic expansion in X/r. The values of t for 
which this expansion breaks down turn out to lie close to 
the limiting value irL 2 , and are thus not very interesting. 
The leading terms of the expansion of v (r) are 

A Tl 

v(r) m — sin (2r + 2(p) + - sin 2 (r + ip). (34) 
2r r 

In this expression 77 is a constant of integration which 
plays a role analogous to the role of C in the one- 
dimensional case - that of a minimization parameter. 
Just as in the one-dimensional case, the minimum of A is 
reached when rj = 0. We will give a perturbative proof of 
this statement at the end of this subsection. With rj = 
the asymptotic expansion of y (r) has the form 



y (r) w cot (r + </?) + — sin 2 (r + ip) . 

4r 

It is convenient to define the parameter 9 as 

e _ ^(0) 
v^(i) 



(35) 



(36) 



in order to cancel out the overall l/y/r dependence of the 
wave function. 9 is always large in the asymptotic region 
t ^> 1. Integrating y (r) in the sense of the principal value 
we obtain 



ln0 = -P 



yMdr + lnOx « l^ln— +ln0i, (37) 
8 r 



where the length scale r is defined as r ~ max (l, |A|), 
and In 9\ represents the contribution to the integral from 
distances r < r . When |A| < 1 this contribution can 
be neglected with logarithmic accuracy. Inverting Eq. 

(p7|) we obtain |A| = 8 In 91 In — . The dimensionless in- 
tegral in the saddle point action J Q v 2 rdr evaluates to 
8 In 2 6 1 In — . and we find 



V (9) ~ exp <^ -Att 2 u 2 D 



In 9 

ln(L/r ) 



(38) 



The envelope of the wave function corresponding to 
the solution described by Eq. ([m]) is 



ip{r) p- (r /r) 



In 0/ In 



(39) 



where we have approximated y/r^^p (ro) ~ ^ (0). For 
ro ~ 1 the error introduced by this approximation leads 
only to O (1) corrections to the large logarithm Int. The 
normalization integral 



N 



27r^ 2 (0) y dr(^) 



2 In 61 In -i- 



(40) 



has two distinct regimes. When 9 2 is greater than 
L/tq the integral is dominated by small distances, and 
the dependence of t on 9 is weak. We will not be 
considering this regime here because it corresponds to 
wave function amplitudes that are close to the limit- 
ing value t = irL 2 . In the opposite case we have 
N = 27r?/; 2 (0)r^L 1 -"/(l-a), where a = hx6 2 /\n^. 

Up to irrelevant constants we then have \ut ~ ln$ 2 , and 
the distribution function becomes 



V (t) ~ exp < —TtviD 



ln 2 t 



In (L/r ) 



(41) 



This result has essentially the same form as the one 
obtained by Falko and Efetov using the tr-model for- 
malism. It differs, however, from their result in two im- 
portant aspects. First, the logarithm in the denominator 
is cut off at distances ~ 1 (pp 1 in conventional units) 
rather than at r ~ I, as in Ref. 0. We believe that this 
difference simply reflects the fact that both the approach 
employed here and the saddle-point solution of the a- 
model represent the same saddle point of the underlying 
theory. However, in order to describe this saddle point 
exactly one has to include short-wavelength degrees of 
freedom which are lost in the derivation leading to the 
tr-model, and that accounts for its inability to obtain the 
correct cut-off scale. 

The second difference may turn out to be an indication 
of a deeper problem with the cr-model. Our calculation is 
performed explicitly for the case when there are no per- 
turbations breaking the symmetry with respect to time 
reversal. Nevertheless, the numerical coefficient in the 
exponent coincides with the answer obtained in Ref. Q] 
for the unitary case (/3 = 2), when the symmetry with 
respect to time reversal is completely broken. In order to 
see what effect violation of time-reversal invariance would 
have in our approach, we explicitly introduce magnetic 
field into the system of equations ( |l5| ) . Choosing the di- 
rection of the field H along the z axis (perpendicular to 
the two-dimensional sample) and writing A = |Hxrwc 
can immediately see that under the assumption of circu- 
lar symmetry of the solution ip, only the terms quadratic 
in A survive. For weak enough fields HL 2 -C <fio (where 
<fro is the magnetic flux quantum) these terms can be 
neglected because they only lead to exponential decay 
of the wave function on a scale l arg er than the sample 
size. Thus the set of equations ( J16| ) is unchanged. On 
the other hand, the fields of this magnitude are sufficient 
to suppress the Coopcron contribution in the cr-model 
calculation leading to a crossover between the orthog- 
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onal and unitary symmetry classes. Therefore our re- 
sults reproduce exactly the cr-model answer in the case of 
broken time-reversal invariance (the unitary symmetry 
class), but they predict a faster decay of the distribution 
function at large t in the orthogonal case. If the fault lies 
with our calculation, it would then mean that while in the 
unitary case a single optimal fluctuation of the potential 
is responsible for the high-amplitude tail of the distri- 
bution of the wave function intensities, an exponentially 
large number of non-circularly-symmetric configurations 
would have to contribute to the probability of finding 
large wave function amplitudes in the orthogonal case. 
(It is certainly very unlikely that accounting for the fluc- 
tuations around the saddle point defined by Eq. (|l^) can 
change the coefficient of a In 2 t term in the exponential) . 
While this possibility cannot be completely discounted 
without a more detailed study of the system Eqs. (15), 
we would like to propose an alternative explanation. 

The expansion into the sum over quasiclassical paths 
that underlies the nonlinear cr-model, and the subsequent 
separation into the diffuson and Cooperon contributions 
implies that self-retracing paths constitute a set of mea- 
sure zero in the overall sum. This is certainly a rigorous 
assumption when there is no "intrinsic" mechanism in the 
calculation that would strongly favor some types of qua- 
siclassical paths over others - as is almost always the case 
in clean chaotic systems where the zero mode treatment 
of the cr-model is applicable. However, when a saddle 
point approximation is being made in order to compute 
the probability of some rarely occurring event, one effec- 
tively selects certain members of the ensemble of random 
potentials that are favorable for such an event. Within 
the cr-model formalism this selection occurrs as a restric- 
tion on the possible configurations of the effective order 
parameter (the Q-matrix, or the supersymmetric density 
matrix q in the formalism of Muzykantskii and Khmelnit- 
skii JT^Q ) , and the information about the actual configu- 
rations of the potential that are being selected is lost. On 
the other hand, as indicated by the calculation presented 
here, these configurations are rather special. Indeed, one 
can interprete the appearance of a quasilocalized state 
possessing an anomalously large amplitude at some point 
r as an opening of a quasi-gap in the spectrum of s-waves 
with respect to r. In a weak periodic potential emer- 
gence of a gap can already be seen in the second order 
of the perturbation theory, and the second order terms 
in the expansion into a sum over quasiclassical paths by 
necessity contain only self-retracing trajectories. Such 
contributions are then overcounted by the cr-model when 
it takes into account the same self-retracing path twice - 
as a diffuson and as a Cooperon. 

While the above argument is rather speculative, it nev- 
ertheless presents a picture which, if confirmed, would 
restrict the applicability of saddle-point solutions of the 
nonlinear cr-model. We believe that this issue deserves 
further consideration. 

To complete the derivation of the distribution func- 
tion presented in this subsection, we outline the proof 



of the statement that 77 = is the optimal choice. The 
solution tp of equations (|I^) can be written as ip (r) = 
^j=- sin (r + (p (r)). Unlike the 77 = case, the phase <p 
does not have a finite limit as r — > 00. The logarithmic 
derivative y (r) of the amplitude function A (r) can be 
expanded in an asymptotic Fourier series of the form 



WW ~EpT j + [y{n,m) cos m (r + p (r)) 

+ y{rL) s ™ m ( r + v(r))}}, (42) 

and ip' (r) can be represented in a similar way. Sub- 
stituting these expansions into Eqs. (|T^) we find that 
2/(i,o) = Thus the choice of the constant of integra- 
tion 77 does not affect the relation between 2/(1,0) and A. 
Integrating y (r) over r we obtain 2/(1,0) hi ■f- & n d there- 
fore the relation between \n8 and A established in Eq. 
([57]) is also ?7-independent to the leading order in A. On 
the other hand, adding the term 2 sin 2 (r + tp) to v can 

only increase the value of the integral Jq v 2 rdr because 
the cross-term in the expansion of the square integrates 
to zero. It follows then that the minimum of the integral 
for a given value of 9 is achieved by setting 77 = 0. Note, 
however, that, as in the one-dimensional case, this proof 
is perturbative - it relics on the possibility to expand the 
solutions in powers of A/r. 

C. The three-dimensional case. 

In the three-dimensional case, it can be self-consist- 
ently demonstrated that large values of t correspond to 
large negative values of A. As a result, there exists a 
range of values of r where expansion in A/r'' -1 = A/r 2 
is impossible. A typical solution which was obtained nu- 
merically using the so-called relaxation method is 
shown in Fig. |IJ One can distinguish three asymptotic 
regimes: (i) r <§; ro, (ii) ro <C r <C r* , and (iii) r ^> r*; 



it will be shown below that ro = 1/|A| and r* = y\X\/2. 
The first region corresponds to a potential well, the sec- 
ond one - to a potential barrier. Taken together, these 
two regions support a resonance in the s-wave channel 
at the energy Ep, However, besides a potential- well - 
potential-barrier combination, resonant scattering can be 
also caused by a weak periodic potential (Bragg reflec- 
tion), and that is exactly what the third region corre- 
sponds to. An interesting consequence of the solution 
presented in Fig. [I] is that the optimal way to achieve 
large values of t in three dimensions is to combine the 
two effects in the "right" proportion. 

Analytically the three regions in Fig. [I] are described 
by the following asymptotic formulae. In the first region, 
v (r) diverges as 



— —2 / 7—2 

tji (r) - A/r + A In r + ( f + c + —A 



(43) 
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where c is an arbitrary constant which cannot be deter- 
mined from the boundary conditions Eq. ((2^). Note that 
the singularity is weak enough so that it produces only 
a finite contribution to the saddle-point action An- The 
behavior of the wave function in this regime is given by 



1 It !t 2 , 

— I — A H — A r mr + cr, 

r 2 3 



(44) 



y(r) 

which corresponds to 

ip(r) &ip{0) (1 + Xr/2 + O (r 2 lnr)) 



The expansion in Eq. ( (43|) breaks down for r ~ 1/|A| 
which gives the approximate value for tq. In the sec- 
ond region the solution can be obtained with the help of 
the quasiclassical approximation, y 2 (r) m v (r) , and the 
result is 



A| * 2/3 



(45) 



Finally, in the third region, where the A/r 2 expansion 
works, v (r) has the asymptotic form 



, . A sin (2r + 2p) sin 2 (r + ip) 
V3 (r) ~ = — — + V \ ^' + O 



(46) 

The constants 77 and 99 in this expression are analogous 
to their counterparts in the two-dimensional case. The 
phase variable tp again has the meaning of the wave func- 
tion phase shift; it is finite due to the rapid decay of the 
potential at large r irrespective of the value of 77. Con- 
stants 77 and p are not independent: they both can be 
regarded as functions of c. Either rj or c can be chosen 
as a minimization parameter. 

We have not been able to establish the analytical de- 
pendence 77 (c); however, based on the numerical analysis 
the qualitative features of this dependence can be de- 
scribed as follows. For a given A there exists a critical 
value Co (A) such that if c > cq, the third (oscillatory) re- 
gion never develops. Instead, v (r) exhibits a singularity 
at some finite value of r, leading to a divergent integral in 
Aq. So the values c > Cq correspond to unphysical solu- 
tions. Exactly at the critical value Co the oscillations are 

2 

also absent, and v (r — > 00) w 1 + A /(4r ). For slightly 
smaller c, 77(c) is large and positive (3> |A|), and the 
oscillations appear only after a more or less protracted 
intermidiate regime in which v (r) ~ 1. An estimate of 
the point r* at which the onset of the oscillatory behavior 
occurrs can be obtained by noting that the oscillations 
of v (r) are driven (through Eq. (|l6|)) by oscillations of 
the wave function. Therefore the amplitude of the os- 
cillations in Eq. (^6|) cannot significantly exceed 1 so as 
to preserve the oscillatory - rather than exponentially 
damped - behavior of the wave function. This require- 
ment leads to 



1 / 2 x 1/4 
r**-(\ +, 2 



(47) 



As c decreases further, 77 monotonically decreases as well, 
eventually covering the whole (+00, —00) interval. 

When 77 is positive or small negative (77 > — |A|) the be- 
havior of v (r) and y (r) in the regions (i) and (ii) depends 
on 77 (or c) only very weakly so that this dependence can 
be ignored in computing the contribution of these re- 
gions to A. Larger negative values of 77 start affecting 
the length of region (ii), and may lead to an emergence 
of a hybrid regime in which v oscillates non-harmonically 
with an r -4 / 3 envelope. 

In what follows we will assume that similarly to the 
two-dimensional case the minimal value Aq is achieved by 
setting rj = 0, even though the breakdown of the asymp- 
totic expansion in — at small distances makes it impos- 
sible to construct an analytical proof of this statement 
analogous to the proofs for one- and two-dimensional 
systems. This assumption is borne out by numerical 
analysis. It is certainly obvious that positive values of 77 
can never be optimal because of a corresponding "costly" 
v ~ 1 region. As for large negative values, they are ruled 
out by the fact that a pronounced "hybrid" regime never 
appears in numerical solutions. Setting 77 = reduces 

Eq. @ to r* » yfi\/2. 

With 77 sa 0, both v-i (r) and 773 (r) reach values ~ 1 
at r = r* . While approximation schemes devised for 
v 1 and v -C 1 break down around r* , we can calculate 



contributions 0\ and 62 to 



_ ^(0) _ 



6162 from regions 



r < r* (i-ii) and r > r* (iii) separately. Similar to the 
two-dimensional case, the L factor in the denominator 
of 9 is introduced to cancel the overall 1/r dependence 
of the wave functions, which is an artefact of spherically 
symmetric boundary conditions. 

(i) - (ii) We combine together the first and the sec- 
ond regions because the wave function amplitude does 
not change appreciably in the very short region (i). Re- 
gion (ii) corresponds to stretched exponential decay of 
•ip. Integrating y (r) = — -\fv~2 (rj over r from 1/|A| to 

A |/2 we obtain 



(iii) Using Eq. ([46]) with 77 = we find 



y (r) ~ cot (r + p) + 



A 

4r 2 



(r + p) 



(49) 



The first term in this expression stems from the oscilla- 
tions of the wave function, while the second one describes 
the decrease (A < 0) of the wave function envelope from 
the center of the sample outwards. The second contribu- 
tion In 62 is given by the principal value of the integral 
/ ydr: 
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In0, 



y (r) dr « r*/4. 



(50) 



Adding the two contributions we establish the relation 
between A and 9: 



4 4 V 2 ' 



(51) 



verifying the self-consistency of the assumption |A| ^> 1 
made at the beginning of this subsection. 

The dimensionless integral in the saddle point action is 
also evaluated separately in the regions (i)-(ii) and (iii): 



f r 192 

/ r 2 v 2 dr — 3r* 3 = — - In 3 9 (regions (i) and(ii)) 
Jo 13 

(52a) 

and 

32 

r 2 v 2 dr = r* 3 /2 = — In 3 6 (region (hi)). (52b) 

Both contributions turn out to be proportional to the 
same power of In 9 indicating that local resonances and 
Bragg reflection play equally important roles in the for- 
mation of anomalously large wave function intensities. 
Combining the results we obtain 



(53) 



The envelopes of the wave function in the two regions 
are 



if} (r) 



V>(0) 



exp{-3r* 2/3 r 1/3 } (for \ < r < r*) 



(54a) 



and 



ij){r) ~ ^^exp{r* 2 /4r- (13/4) r*} (for r > r*). 



The normalization integral is given by 

1 



N = (0) \ ^ + 4£e- 



(54b) 



(55) 



When L > 9 2 (with more realistic boundary conditions 
this inequality will become L 3 ^> 9 2 ), the normalization 
integral is dominated by the contribution from region 
(iii). Then ip 2 (0) oc 9 2 , and we finally obtain with expo- 
nential accuracy 



V (t) ~ exp 



7 



2197 



(p F l)\n 3 t 



(56) 



Of course, the separation into regions r < r* and r > r* 
is approximate. It is possible that the crossover region 
r ~ r* gives a contribution of the same order of magni- 
tude as the two asymptotic regions (ii) and (iii). Thus 
the number 7/2197 « 3.2 x 10~ 3 can only be consid- 
ered as an order-of-magnitude estimate of the coeficicnt 
k introduced in Section |J. It must be mentioned that 
convergence to the asymptotic form of Eq. ( |56| ) is ex- 
tremely slow because of the stringency of the requirement 
r*=j|ln0»l. 

In the opposite regime L <C 9 2 , tp 2 (0) is close to its 
maximal value ~ 1, and therefore its 9 dependence is 
very weak. On the other hand, similarly to the one- 
dimensional case, 9 2 becomes proportional to the electric 
response time t, leading to the distribution of t having 
the form identical to Eq. (|56|), 



V (t) ~ exp {-k (p F l) In tA} 



(57) 



Both Eq. ([36() and Eq. (|5^) differ significantly from 
the corresponding cr-model results 0,^]. We discuss the 
possible origins of this difference in the next section. 



IV. DISCUSSION 

A. General remarks. 

The main result of the work presented in this paper is 
that statistics of seldom occurring events in disordered 
electronic systems can be successfully studied using a 
variant of the optimal fluctuation method. All previous 
approaches to the problem have been based on various 
formulations of the nonlinear cr-model, and they invari- 
ably seem to require an extension of the cr-model to, and 
sometimes beyond, its limits of validity. In retrospect, 
this situation seems rather natural. Indeed, the success of 
the cr-model in describing a wide variety of phenomena in 
chaotic and disordered systems can be traced to the fact 
that most such phenomena are quasiclassic in nature and 
therefore their most appropriate description is in terms 
of an expansion into sums over quasiclassical trajectories. 
Calculation schemes based on non-linear cr-models allow 
one to perform such summations in a very efficient man- 
ner by introducing collective variables (Q-matrices) that 
change slowly in space [|],[|. When low-order moments 
and correlation functions are described the integrals over 
Q-matrices get almost equal-weight contributions from 
the whole degenerate (or near-degenerate) manifold 
of allowed Q-matrix configurations. This effectively cor- 
responds to probing a large number of disorder configu- 
rations. 

An investigation of the statistics of rare events, on the 
other hand, presents quite a different type of problem. 
The configurations of the random potential that give rise 
to such events come from a small subset of all possible 
configurations. This circumstance was reflected in the 
approach developed in Ref. M - the large parameter t 
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served to strongly break the degeneracy of the realiza- 
tions of the Q-matrices, and only the vicinity of a saddle 
point of the integral over Q's was shown to give a signifi- 
cant contribution to the probability of observing large t. 
What was effectively encoded in such a scheme was the 
fact that the purpose of the calculation was to select those 
special potentials that are capable of producing a given 
value of t, and a corresponding selection of the relevant 
configurations of Q's was a computational tool allowing 
to do that. 

The disadvantages of such an approach are obvious: 

(i) First, the assumption that the motion of elec- 
trons can be described entirely in quasiclassical terms 
imposes certain restrictions on the types of the poten- 
tials over which averaging is performed. As follows from 
the calculations presented here, in three spatial dimen- 
sions there exists a region r ~ r* where the solution of 
the Schrodinger equation cannot be obtained in the qua- 
siclassical approximation. As a result, the tr-model fails 
to recognize the existence of this region and it also misses 
entirely the local resonance formed at r < r* . 

(ii) Second, even when the er-model calculations suc- 
ceed in correctly - albeit implicitly - identifying the 
relevant disorder configurations (as in one- and two- 
dimensional samples) they do not always produce correct 
answers because contributions from the short- wavelength 
degrees of freedom (massive modes) eliminated in the 
transformation from the fast to slow variables are missed. 



B. The direct optimal fluctuation method. 

It is not immediately evident, however, that a rela- 
tively naive approach based on the direct search for an 
optimal fluctuation should be more reliable. In order for 
it to work, the probability of observing a large value of 
t must be determined by a sum over disorder configu- 
rations which all come from a single compact region of 
the configurations space. The potentials forming this re- 
gion differ only slightly from some optimal configuration 
which corresponds to the saddle point. Although we have 
not proved in this work that the saddle point identified 
here gives the dominant contribution to the functional 
integral, "the preponderance of evidence" based on the 
comparison of our results with those obtained using the 
cr-model would indicate that this is indeed the case. 

The starting point of our qualitative analysis is the 
observation that, apart from the difference in the cutoff 
scale, our variant of the optimal fluctuation method (the 
"direct" optimal fluctuation method) reproduces identi- 
cally in the two-dimensional and one-dimensional cases 
the tr-model results for the unitary ensemble. Assum- 
ing that more than a chance coincidence is involved, it 
is reasonable to conclude that the optimal configurations 
of disorder found in this work are the same as the ones 
that are responsible for the saddle point of the cr-model. 
In other words, there must exist a one-to-one correspon- 



dence between the saddle points of the theory defined by 
Eqs. ( |l4| ) and the saddle points of the cr-model. Assum- 
ing the existence of such a correspondence we will try to 
elucidate the origins of the cr-model results by analyzing 
the corresponding optimal configurations of the potential 

Before proceeding with this analysis we would like 
briefly to address the question of the stability of the 
saddle point described by Eqs. (fl5|). In the three- 
dimensional case the dominant contribution to the saddle 
point action Ao comes from small distances where the 
optimal potential is much larger than its typical Gaus- 
sian fluctuation, and stability with respect to fluctuations 
does not pose a problem. In two dimensions, however, 
the outlying regions of the optimal configuration - where 
the magnitude of the potential vanishes as 1/r — must be 
taken into account. In order to argue that fluctuations 
do not destroy the saddle point we note that although 
we find it convenient to write the Gaussian distribution 
function W [U] for the potential in the coordinate rep- 
resentation, it can be written in any orthonormal basis 
{fn ( r )}- A typical amplitude of a dimensionless basis 
function /„ in a typical configuration U (r) is l/y/v^D. 
Let us now choose one of the basis functions, say fo (r), to 
be proportional to the optimal solution v (r) given by Eq. 
(p4|). Using the normalization condition J drf^ (r) = 1 
we obtain 



and 



/o(r) 



v (r) 



1 



sin 2r 



7rln(L/ro) T 
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Vln(L/r ) 



/o(r) 



(58a) 



(58b) 



Thus the optimal fluctuation has a much larger ampli- 



tude than the typical one as long as In 9 



l n(L/r ) 

This condition, of course, is just a natural requirement for 
the validity of the optimal fluctuation method | In V | ^> 1 . 

It remains to be shown, however, that other compo- 
nents of a typical fluctuation U (r) (i.e. those orthogonal 
to fo) do not destroy the saddle point. A rigorous in- 
vestigation of the fluctuations around the saddle point 
defined by Eqs. (|l5|) will be the subject of a forthcoming 
publication. Nevertheless, a plausible argument in favor 
of the stability of this saddle point can be made based 
on the following observation. Since the appearance of 
anomalously localized states due to Bragg reflection is a 
phenomenon analogous to the emergence of a band struc- 
ture in a periodic lattice, the suppression of this effect by 
fluctuations of the potential is equivalent to the local- 
ization transition which destroys the band structure in 
ordinary periodic lattices. Therefore L < L c seems to be 
a sufficient condition for the stability of the saddle-point 
solution in the two-dimensional case. 
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C. The asymptotics of the distribution functions. 

1. The one- dimensional case. 

We will now try to use the physical intuition afforded 
by the optimal fluctuation concept to compare the cr- 
model results of Eqs. (0) and (|^) with the distribution 
functions Eqs. (f|), @, @ and @ derived in the 
preceding Section. The first question that can be easily 
answered is why the distribution function for the electric 
response times in the one-dimensional case is log-normal 
instead of a power law (i.e., exp (— C\ lni)) obtained by 

a naive extrapolation of the exp d ln d tj dependence 

to d = 1. In one and two dimensions the optimal config- 
urations found in the preceding Section are "global" , i.e. 
the integration in the saddle-point action Jv 2 dr must 
be extended over the whole sample. In contrast, in the 
three-dimensional case the optimal fluctuation, even with 
the oscillating tail included, is local, so that the above in- 
tegral converges at large distances. It is well known |21[| 



that distribution function tails of the e 



■ In" t 



type usu- 



ally appear as probabilities of optimal fluctuations con- 
fined to a finite volume. Thus in order to explain the 
"anomaly" in the one-dimensional case we have to un- 
derstand why a local fluctuation of the random potential 
necessary to achieve a given value of t has a lowerprob- 
ability than the " gl obal" one proposed in Ref. || and 
rederived in Section III A . A local fluctuation of the po- 



tential leading to a large value of t would have to be able 
to support a narrow resonance. Assuming as in Section 
[n| a rectangular shape for the potential barrier responsi- 
ble for the formation of the resonance we can repeat the 
calculation presented there almost verbatim except that 

t must everywhere be replaced by tA. We then obtain 
2 

In? k — — (ppl)\ntA. In order for this behavior to 
3v3 

dominate we must have ppllntA <C (//L)ln 2 fA, or 

lniA>p F L. (59) 

However, from Eq. (^) we see that the corresponding val- 
ues of b - which determine the size of this local fluctua- 
tion of the potential - become larger than the length of 
the sample 2L. This is clearly unphysical. Thus despite 
a slower i dependence of the probability of local reso- 
nances, the same value of t has a much larger probability 
to be produced by an accidentally formed Bragg mirror 
for all reasonable values of t. 

Note also that IniA ~ ppL corresponds to v ~ 1 (or 
U ~ Ep), which invalidates the perturbative expansion 



of Section III A Essentially, the values of IniA of the 
order of ppL or larger correspond to a trivial case: a 
sample is insulating because the potential is larger than 
Ep almost everywhere except for a small island in the 
center where almost all the weight of an eigenstate at 
Ep is concentrated. In this regime the distinction be- 



tween "global" and local fluctuations of the potential gets 
blurred. 

It is interesting to note that a calculation based on 
the ballistic cr-model performed in Ref. || for quasi- 
one-dimensional conductors indicates the existence of a 
crossover from log-normal to power-law distribution at 

In t A ~ — . Although we have not investigated the quasi- 
one-dimensional case here, it is likely that the physical 
picture of the interplay between the "global" and local 
fluctuations discussed above should not be much differ- 
ent. If that is the case, then it is probable that the 
crossover found in Ref. H has as its underlying cause 
the same mechanism of local fluctuations becoming com- 
parable in size to the length of the sample. This hypoth- 
esis, however, leaves unexplained the difference in scales 

at which the crossover occurrs — - in Ref. j|] as opposed 

to ppL in the argument presented above. It is possi- 
ble that the quasi-one-dimcnsional case brings in some 
new features that are not recognized by the estimates 
based on the purely one-dimensional model. It should be 
mentioned, however, that none of the variants of the cr- 
model can provide an adequate description of the effects 
associated with local resonances, and it is possible that 
a cr-model estimation of the crossover scale may not be 
entirely reliable. 



2. Two-dimensional conductors. 

A rather comprehensive discussion of the results ob- 
tained by the direct optimal fluctuation method in the 
two-dimensional case and their counterparts established 
using the cr-model formalism has already been presented 
in Section [II B . Here it seems appropriate to briefly re- 



iterate the following two main points of that discussion. 
First, it comes as no surprise that when massive modes 
are taken into account, the short distance cut-off scale in 
the logarithm determining the system size dependence of 
h\V becomes of the order of the electron wavelength p'p 1 
rather than the mean free path I. This leads, however, 
only to a small relative change in the cr-model result. 

On the other hand, the apparent ensemble indepen- 
dence of the distribution function asymptotes obtained 
by our method stands in a striking contradiction with 
the letter and the spirit of the cr-model and is quite enig- 
matic. 



3. The three-dimensional case. 

The same as in the one-dimensional case issue of 
whether the cr-model calculations can correctly handle 
the role of fluctuations of the potential with the ampli- 
tude larger than the Fermi energy arises in connection 
with the In 2 t asymptotics of the distribution of wave 
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function amplitudes obtained in Ref. for the three- 
dimensional case. Derivation of the cr-model involves lin- 
earization of the spectrum near the Fermi energy. As 
a result, there is no intrinsic scale in the model that 
would relate the amplitude of the fluctuating potential 
to the electron energy. The only scale is provided by 



1 



the dispersion of the fluctuations of the potential 

TTV d T 

This limitation of the model does not affect the computa- 
tion of probabilities of typical events, or of the averages 
that are dominated by such events, because typical po- 
tentials are small compared to the Fermi energy. How- 
ever, a problem arises when rare large-amplitude fluc- 
tuations of the potential U (r) become dominant. The 
cr-model cannot detect the existence of classical turn- 
ing points around which the quasiclassical approximation 
breaks down. Thus a possible explanation of the result 
of Falko and Efetov M is that within the cr-model ap- 
proach the classical turning point at r* is missed, and 
the Bragg mirror is effectively assumed to persist un- 
til distances of the order of the mean free path Z. In 
the notation of Section III C| that would correspond to 

— IAI 2 

|A| ~ llnd, and then Aq oc J t v 2 r 2 dr oc cx / In 2 6, 

leading to InV ~ — (ppl) 2 In 2 t which coincides with the 
answer obtained by Falko and Efetov. The reason that 
such a cut-off scheme leads to a higher estimate for the 
probability V is that it corresponds to an incorrect solu- 
tion of the Schrodinger equation in the region of large po- 
tentials U > Ep which overestimates the rate of growth 
of the wave function amplitude towards the center of the 
sample. 

Another discrepancy between our results and those ob- 
tained with the help of nonlinear cr-models in the three- 
dimensional case is the difference in powers of (pfI) in 
the exponents in Eqs. (Q) and (|57|). We believe that 
this discrepancy has the same origin as the difference 

between In — and In — in the two-dimensional case - 
I r 

the error introduced by using the mean free path Z to 
determine the cut-off scale. It is interesting to note that 
while the cut-off procedure used in the diffusive cr-model 
approach of Ref. j|) is capable of producing only an 
order-of-magnitude estimate C3 ~ (pfI) 2 , the ballistic 
cr-model of Muzykantskii and Khmelnitskii |l(]||, while 
giving an illusion of computing the coefficient C3 exactly 

(pfO )) nevertheless leads to the same extra 



(C 3 



9\/3 



power of (pfI). To explain this seemingly paradoxical 
situation we will first examine the cut-off procedure em- 
ployed in Ref. ||. It is based on the condition, pointed 
out in Ref. Q , that in order for the diffusive a- model to 
be applicable, the spatial gradients of the Q-matrix com- 
ponents cannot exceed 1/Z. Explicitly, for the calculation 
performed in Ref. H this condition reads 



where Ai parametrizes the non-compact bosonic sector of 
the Q-matrix || (not to be confused with the Lagrange 
multiplier A used throughout this work). The distance 
Z» where this condition is violated is used in Ref. || as 
a short-distance cutoff. It was conjectured in Ref. R] 
and later confirmed in Ref. (20) that spatial structure 
of the saddle-point solution for Ai mimics the envelope 
of anomalously localized states described by the saddle 
point of the cr-model. Therefore the optimal configura- 
tion of (d/dr) In Ai (r) corresponds to the non-oscillating 
part of y (r) in the direct optimal fluctuation method. 

This correspondence allows us to see directly what ef- 
fect an artificial short-distance cutoff at Z* would have 
in our approach. From Eq. ( |49| ) we find Z? ~ Z|A|. In- 
troducing such a cut-off at Z» into Eq. (|50| ) we obtain 
|A| ~ Z* hi# and therefore Z* is estimated as Z* ~ lh\0. If 
then the integral determining the saddle-point action Aq 
is also cut off at Z* , we obtain 



.4,, - / / v 2 r 2 dr ~ Z^ 



dr 



lnAi 



< 



1 

V 



(60) 



Z 2 In 3 0, 



or (pfI) 2 m3 ^ m conventional units, leading to the incor- 
rect answer that was already quoted in the Introduction. 

It is thus evident that distances shorter than Z* give 
an important contribution to the saddle-point action, and 
this contribution cannot be accounted for by the diffusive 
nonlinear cr-model. It should be emphasized that exclud- 
ing the short-distance contribution in this way leads to 
a significantly smaller estimate for the probability of ob- 
serving a given (large) value oft. This can be understood 
by using the following argument. The wave function am- 
plitude grows substantially between r ~ Z* and r ~ 1. 
Neglecting this growth leads to a need for a faster in- 
crease in t between L and Z* which can only be achieved 
by means of a "costly" boost in the amplitude of the 
Bragg mirror. As a result, non-optimal configurations of 
the random potential are selected. 

Turning now to the calculation performed with the 
help of ballistic nonlinear cr-model for the three-dimen- 
sional case in Ref. || we notice that the distance r* which 
separates the "reaction" and "run-out" zones is of the 
same order as Z*. The contribution to the saddle point 
action from the "run-out" (diffusive) zone leads to the 

7T 2 

already quoted ~^j= (pfI) value for C3, while the "re- 
action" zone produces a contribution that has one less 
power of the large logarithm In 9 and is thus negle cted. 
In contrast, in the calculation presented in Section III C 
of this paper, the contribution from distances of the or- 
der of r* ~ \n9 <C Z* dominates the saddle point action 
and leads to a larger estimate (Eq. (|5a|)) for the proba- 
bility of observing anomalously high values of t in the 3- 
dimensional case than the one obtained in Ref. Thus, 
the ballistic generalization of the nonlinear cr-model is 
also not capable of detecting the existence of the scale r* 
at which the quasiclassical approximation breaks down. 
Moreover, since the calculation within the framework of 
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the ballistic cr-model does not involve any ultraviolet di- 
vergencies that would necessitate a short-distance cut-off 
as in the case of the diffusive cr-model, it appears plausi- 
ble that the ballistic variant of the model is equivalent to 
introducing an ultraviolet regularization into the theory. 



D. Prelocalized wave functions. 

To complete the comparison of the results obtained by 
the direct optimal fluctuation method and those found 
using the nonlinear cr-models we now turn to the issue of 
the shape of the envelope of anomalously localized states 
that are responsible for the large-i tails of the distribu- 
tion functions V (t). In Ref. f?J it was found that in two 
dimensions such states are characterized by amplitudes 
decaying outwards according to a power law: 



lnt/ln(L/Z) 



(61) 



This result has to be compared with Eq. which can 
be rewritten as 



hKr)r 



V> 2 (o) /r \ lnt/ln t 



(62) 



Apart from the difference in the cut-off scale (I vs. ro) 
which was discussed at length above, Eq. (|6^) contains 
an extra 1/r factor in the denominator. It is simply a 
consequence of the idealized model adopted here with its 
circularly symmetric boundary conditions. In a more re- 
alistic model the optimal wave function would become a 
superposition of different angular momentum eigenstates. 
The l/v^ behavior of the circularly symmetric compo- 
nent would be cancelled in such a superposition. Taking 
into account fluctuations around the saddle point would 
also have the effect of suppressing the (1/^/r) 2 factor in 
the spatial dependence of the averaged envelope of the 
optimal solution. Note also that Eq. ( |6l| ) does indeed 
describe the averaged envelope of anomalously localized 
states. 

In three-dimensional samples the states with anoma- 
lously high local amplitudes were found in Ref. to have 
the following envelope: 



|^(r)r~e X p<j-A(l-- 



(63) 



where A is a constant which in the leading logarithmic 
approximation is equal to hit. Comparing this to Eq. 
(p4) we see, in accordance with the discussion above, 
that the cr-model result gives the correct functional form 
e Const./r on jy f or t ne region of large r which corresponds 
to the region of space where the optimal potential forms 
a Bragg mirror. The estimate ilnt for the constant in 
the exponential, however, contains an extra factor of I 
compared to Eq. (54b) which is a consequence of the 



choice of the cut-off scale made in Ref. |7j]. The region 
of the stretched exponenti al de cay of the wave function 
envelope described by Eq. ( 54a ) is missed in the cr-model 
calculation entirely. As in the two-dimensional case, the 
extra powers of r in the denominators of Eqs. (|w|) are a 
consequence of the artificial rotational symmetry. 

It was conjectured in Ref. j?J that these high-amplitude 
states have a complicated "snake-like" structure at short 
distances. The conjecture was based on the fact that the 
method employed in Ref. was applicable to amplitudes 
as high as 



t< 



Vp d 



(64) 



V 



rather than a naive t < expected from a cut-off at I. 

We have found no evidence for such behavior here. Since 
the solutions of the saddle point equations were assumed 
to be rotationally invariant from the outset, this cannot 
be regarded as a conclusive evidence against such a sce- 
nario. However, potentially more important is the fact 
that we did not encounter any limitations on the possible 
values of t analogous to Eq. (|64|) . Analysis of the fluctua- 
tions around the saddle point defined by Eqs. ( |l5| ) as well 
as an investigation of the possibility for non-rotationally- 
invariant solutions of these equations is needed to settle 
the issue conclusively. 



E. Are the cr-model results universal? 

Finally, we would like to make a few remarks concern- 
ing the issue of universality. An important consequence 
of the dominant role played in three dimensions by large 
local fluctuations of the potential that are responsible 
for the formation of local resonances is the non-universal 
character of the distribution functions derived in Section 



IIIC. Indeed, such large-amplitude configurations of the 
potential can only be optimal if they are not too "expen- 
sive" - i.e. if their action A is not too high - compared 
with the low-amplitude "global" fluctuations of the po- 
tential (Bragg mirrors). Generalizing the distribution of 
the random potentials to 



W n [U] oc exp 



I' 



(r) d d r 



(65) 



where a is the dispersion and n > 1, we find through a 
model calculation similar to the one performed in Section 
H that forming states with large amplitudes at some r by 
means of local resonances always leads to \riV oc — In t. 
For n < d/2 the "cost" of a corresponding Bragg mir- 
ror can be estimated from an appropriate generalization 
of Eqs. (|l~6|), and it also leads to a \n d t dependence, 
as was demonstrated by a detailed study of Gaussian 
(ri = 1) distribution in three dimensions. Thus local 
large-amplitude fluctuations of the potential are always 
important for small enough (< d/2) values of n. In the 
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marginal case n = d/2 the weight of the Bragg mirror is 
given by 



InV 



ln d t 



In- 



L 

ro 



d-l ' 



(66) 



making it a much more probable way of achieving a large 
local wave function intensity t. When n > d/2 the prob- 
ability of the corresponding Bragg mirror has a faster 
than ln d t dependence on lni, but it is compensated by a 
power of the system size in the denominator: 



ln'P cx 



ln 2 "t 

Jjln—d ' 



(67) 



This situation is realized in the Gaussian case when d = 
1, and it was already discussed in detail. Therefore the 
seemingly universal character of the log-normal behavior 
of the tails of the distribution functions is tied to the 
assumption - which is a starting point in the derivation of 
the er-model - that the random potential has a Gaussian 
distribution. 



V. CONCLUSIONS AND OPEN QUESTIONS. 

The most important problem arising from the results 
of the present study is the need to explain the fact that 
the asymptotics of the distribution functions derived us- 
ing the direct optimal fluctuation method do not exhibit 
any dependence on weak magnetic fields. This feature 
of our answers must be contrasted with all the previous 
calculations performed in the framework of nonlinear a- 
models in which it is quite obvious that when Cooper 
modes acquire a mass as a consequence of broken invari- 
ance with respect to time reversal, the number of inde- 
pendent components of the Q-matrices changes and that 
has a profound impact on the results. There are two 
alternative resolutions of this contradiction. 

One possibility is that our assumption that rotationally 
invariant solutions of the saddle-point equations domi- 
nate the saddle-point action is not valid in the absence 
of magnetic field. A more detailed study of the properties 
of Eqs. (|l5|) will be the subject of a future publication. 

The second possible explanation is that there is some 
internal inconsistency in the cr-model which causes it in 
certain circumstances to overcount the number of degrees 
of freedom. In our view, this possibility must be taken 
seriously. 

Whether or not the fluctuations around the saddle 
point can change the leading order terms in mi is also 
one of the questions that were left beyond the scope of 
this work. A peculiar disagreement in the crossover scale 
from In 2 i to In i asymptotics of In V between the purely 
one-dimensional and quasi-o nc-dimcnsional cases which 
was noted in Section IV C 1 makes it rather desirable to 



extend our method to the quasi-lD and quasi-2D geome- 
tries. 

To conclude, the main results of this work can be sum- 
marized as follows. We have demonstrated that the opti- 
mal fluctuation method is a useful tool for the investiga- 
tion of the statistical properties of anomalous electronic 
eigenstates in disordered two- and three-dimensional con- 
ductors. It was shown that in three dimensions this 
method is preferable to the nonlinear cr-model because 
the latter does not include effects associated with local 
resonances that can be formed by the random potential. 
In the two-dimensional case the results obtained by the 
optimal fluctuation method essentially coincide with the 
cr-model results for the unitary ensemble of random po- 
tentials which we interprete as an indication that the sad- 
dlepoint of the reduced nonlinear cr-model found in Refs. 
lEim!] corresponds to the same saddle point of the full 
theory as the one that describes the optimal fluctuation 
of the potential. 
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APPENDIX: 

Introducing an auxiliary field \ as m the main text, we 
can represent M^, as 



j 'v^V (£) e 'J ^0 



T +U(r)-E U« 



(Al) 



Performing a 7r/4 rotation in the (x>"0) space 

X=(V^i+^2)/V2, 4>=(il>i-ih)/y/2, (A2) 



and introducing infinitesimal convergence factors we ob- 
tain 



V 



Vi 



v 



V'2 



x e 



2wJ Vv2ti\ 

§i J drfViO'+'Vi-^O'-'iM 



(A3) 



where — §^+U (r) — E±iO. Thus, up to a constant, 
Afip is equal to a symmetrized spectral determinant: 



Af,/, cx det | E - H + io\ det | 

_ e 3CTrln(E-j?+tO) 



E-H-iO\ 



(A4) 



The variation of the logarithm in the exponent with re- 
spect to U (r) is equal to the real part of the Green's func- 
tion at coinciding arguments 5iG (r, r; U) . This quan- 
tity probes the whole band and is not sensitive to small 
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changes in U. It therefore can be treated as a constant 
(which we denote as G), so that M^, becomes 



K r G f drllfr 

Alip oc e J 



(A5) 



The background potential J drU (r) can be absorbed into 
a redefinition of energies, which justifies the approxima- 
tion made in deriving Eqs. (plf). 

It should also be noted that in the sigma model for- 
malism the real part of the Green's function G (r, r) is 
effectively set to zero under the assumption of an infinite 
symmetric band. Therefore any corrections to the distri- 
bution function that may arise due to Af^ ^ 1 are at any 
rate beyond the scope of the cr-model. 



[1] Mesoscopic Phenomena in Solids, B. L. Altshuler, P. A. 

Lee and R. A. Webb, eds., North Holland 1991. 
[2] B. L. Altshuler, V. E. Kravtsov and I. V. Lerner, in 0, 

and references therein. 
[3] F. Wegner, Z. Phys. B 35, 207 (1979). 
[4] B. A. Muzykantskii and D. E. Khmelnitskii, Phys. Rev. 

B 51, 5480 (1995). 
[5] K. B. Efetov, Adv. Phys. 32, 53 (1983). 
[6] V. I. Falko and K. B. Efetov, Europhys. Lett. 32, 627 

(1995). 

[7] V. I. Falko and K. B. Efetov, Phys. Rev. B52, 17413 
(1995). 

[8] A. D. Mirlin, Phys. Rev. B53, 1186 (1996). 

[9] B. A. Mu zykantskii and D. E . Khmelnitskii, unpublished 
(preprint cond-mat/9601045). 
[10] B. A. Muzykantskii and D. E. Khmelnitskii, JETP Let- 
ters 62, 76 (1995), [Pis'ma Zh. Eksp. Teor. Fiz. vol. 62, 
p. 68 (1995)]. 

[11] B. Altshuler and V. Prigodin, JETP Letters 47, 43 
(1988), [Pis'ma Zh. Eksp. Teor. Fiz. vol. 47, p.36 (1988)]. 

[12] B. I. Halperin and Melvin Lax, Phys. Rev. 148, 722 
(1966). 

[13] J. Zittartz and J. S. Langer, Phys. Rev. 148, 741 (1966). 

[14] Such a formulation of the problem implies perfect ro- 
tational symmetry and sidesteps the issue of the role 
of irregular boundary conditions. However, as will be- 
come evident from the answers, the asymptotic form of 
the three-dimensional distribution is not sensitive to the 
boundary conditions at all. The two-dimensional case is 
trickier and may require a separate study of the stability 
of the solution with respect to the change in boundary 
conditions. Nevertheless, since the a-model calculations 
were performed for similar geometries, the model adopted 
here is completely adequate for the purposes of compar- 
ing different approaches. 

[15] In the formalism of Falko and Efetov J(| the extra delta- 
function in the definition of V was essential in enabling 
the mapping of the problem onto the nonlinear <r-model. 
However, the saddle-point solution of the reduced variant 
of the model which was derived in Ref. || to describe de- 



[16] 
[17] 

[18] 

[19] 

[20] 
[21] 



viations from the universal Porter-Thomas statistics |l(| 
was shown to involve only the variables parametrizing 
the bosonic sector of the theory. This can be viewed as 
an indication that effects associated with the statistics of 
energy levels are not important in the approximation cor- 
responding to the saddle-point treatment of the nonlinear 
a-model. In the context of the direct optimal fluctuation 
method this observation receives quite a natural expla- 
nation. A given value of energy E can be fine-tuned to 
become an eigenvalue by a shift in the optimal config- 
uration of U (r) of the order of the mean level spacing 

A = — — . Such a shift can result in corrections to InV 

that are at most ~ Eft - much less then the leading 
contribution when kit is large. 

Haake, Quantum signatures of Chaos, Springer- Verlag, 
Berlin, New York 1991. 

In several cases the probabilities obtained by our method 
are smaller than the correponding a-model results. The 
relation between the two sets of answers and the likely 
reasons for the differences between them are discussed in 
Sections [in] and |jv[ 

E. Kamke, Differential Equations, Chelsea Publ., NY 
1971. 

W. H. Press et.al, Numerical Recipes in C, Cambridge 



University Press, Cambridge 1988. 

A. D. Mirlin, unpublished (preprint cond-mat/9512095) 
I. M. Lifshitz, Soviet Phys. - JETP 17, 1159 (1963). 



FIG. 1. The profile of an optimal configuration of the po- 
tential in three dimensions as a function of the radial coordi- 
nate. 



16 



0.4 




p F r* 5 10 15 20 

P F r 



Fig. 1 



